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Abstract 



The mass distributions of baryon resonances populated in near-central 
collisions of Au on Au and Ni on Ni are deduced by defolding the p t spectra 
of charged pions by a method which does not depend on a specific resonance 
shape. In addition the mass distributions of resonances are obtained from 
the invariant masses of (p, n^) pairs. With both methods the deduced mass 
distributions are shifted by an average value of —60 MeV/c 2 relative to the 
mass distribution of the free A(1232) resonance, the distributions descent 
almost exponentially towards mass values of 2000 MeV/c 2 . The observed 
differences between (p, tt~) and (p, tt + ) pairs indicate a contribution of 
isospin I = 1/2 resonances. The attempt to consistently describe the 
deduced mass distributions and the reconstructed kinetic energy spectra 
of the resonances leads to new insights about the freeze out conditions, i.e. 
to rather low temperatures and large expansion velocities. 



1 Introduction 

In central heavy ion collisions the number of mesons produced increases with 
bombarding energy. In the energy range between 1 and 2 AGeV, 8% to 22% of 
the available energy is dissipated into this channel which predominantly consists 
of pions |], |2|. The amount of dissipated energy also depends on system mass: 
In central Au + Au reactions at 1 AGeV it is 8%, but 15% in central Ni + Ni 
reactions at the same bombarding energy. The process responsible for meson 
production is believed to be predominantly the excitation of baryon resonances 
during the early compression phase of the collision. In the later expansion phase 
these resonances decay. The primordial decay yields a number of different mesons, 
mainly pions. Mesons with masses higher than the pion mass, e.g. 7], p, may 
subsequently decay also into several pions. In total, the average mass of the 
excited baryon resonances and the number of pions produced by their decay 
chains increase with bombarding energy. This mechanism is the basic process 
of pion production which is used in nuclear transport models to describe the 
dynamics of relativistic heavy ion collisions 0, f|, [|, ||, [7], ||, [9|, [R], |TT| . In a recent 
publication [|ll[] it was shown that all resonances up to the A(1950) resonance have 



to be included to successfully reproduce the experimental pion data at energies 
between 1 and 2 AGeV. 

The experimental confirmation that baryon resonances are excited in rela- 
tivistic heavy ion collisions is rather scarce. The earliest indication came from 



the comparison by Brockmann et al. [12| of the pion and proton energy spectra 



observed in 1.8 AGeV Ar + KC1 reactions. This comparison was later extended 
with the attempt to determine the A(1232) resonance mass and temperature 



131. Similar analyses of the pion energy spectra measured in 1 AGeV Au + Au 



collisions were published in [14|. At much higher energies, i.e. 13.7 AGeV, the 



fraction of nucleons excited to the A(1232) resonance was determined for Si + 



2 



Al, Pb reactions fI5 |. In this study also the invariant mass distribution of the 



A(1232) resonance was reconstructed from correlated (p, n + ) pairs. 

To identify structures in the invariant mass distribution of correlated proton 
and pion pairs provides the direct proof that nucleons are excited to high-lying 
resonances. The major obstacle in reconstructing the invariant mass is the large 
background of non-correlated (p, 7r) pairs, which increases with the number of 
protons and is therefore particularly large for central reactions of heavy-mass 
systems. In peripheral reactions with very light projectiles, e.g. p |T6| , |17jj or 
3 He |18| induced reactions at around 2 GeV bombarding energy, the (p, tt + ) cor- 
relations were successfully analyzed and the mass distribution of the A(1232) 
resonance was determined. The resonance mass was found to be shifted by —25 
MeV/c 2 to lower masses in reactions on various targets, compared to those on 
protons. Very recently the A (1232) was reconstructed from (p, tc + ) pairs detected 
at 1.97 AGeV bombarding energy in 58 Ni + Cu collisions |TjJ. In this case much 
larger mass shifts of —75 MeV/c 2 were observed in central collisions, the shift 
became smaller with increasing impact parameter. That the (p, tt ) pairs can 
also be successfully studied in central collisions of very heavy systems at 1 AGeV 



bombarding energy was shown in ||20|| . This report contained a preliminary anal- 
ysis of the data used in the present publication. Although at 1 AGeV energy 
the excitation of the A(1232) resonance is the dominating channel there is direct 
experimental evidence from the identification of the i] meson [BT], that at this 
energy also higher-lying resonances, in particular the iV(1535) resonance, must 
become excited. If this is the case then the resonance decay into the 2tt channel 
also has to be present, and the analysis of (p, 7r ± ) correlations will not yield a 
complete picture of the resonance distribution. On the other hand, the 2tt decay 



channel can be observed in the pion transverse momentum spectra fTT 



The subsequent sections are organized as follows: The section 2 contains a 
short description of the experimental procedures. A more detailed description 
can be found in The methods to extract the invariant mass distributions 
of the resonances involved are presented in section 3. These methods include 
the analyses of the measured ir~ and 7r + transverse momentum spectra, and 
the analyses of the (p, tt~) and (p, tt + ) correlations. In both cases the residual 
Coulomb interaction between the charged pions and the positively charged baryon 
distribution causes a perturbation fl], ||, which is of opposite sign for 7r~ and n + 
and which we compensate by analyzing the average of the ir~ and 7r + spectra, 
respectively by taking the average invariant mass distribution from (p, tt~) and 
(p, 7r + ). The justification to apply this procedure is given in the appendix 6.1. The 
experimental results are presented in section 4, section 5 contains their discussion 
and the summary. 
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2 Experimental procedures 



The reactions 197 Au + 197 Au and 58 Ni + 58 Ni were studied at nominal beam 
energies of 1.06, respectively 1.06, 1.45 and 1.93 AGeV. The 197 Au and 58 Ni 
beams were accelerated by the UNI LAC / 'SIS accelerator combination of the 
GSI /Darmstadt. The duty cycle was 75% with a spill length of 4 s. The av- 
erage beam intensities varied between 1 • 10 5 and 5 • 10 5 particles per spill. The 
targets consisted of self supporting foils with a thickness of 100/xm(Au), respec- 
tively 270/xm(Ni). This corresponds to an interaction probability of 0.5%. The 
energy loss of the beams in the targets is of around 0.01 AGeV and was neglected. 

The particles produced by the Au + Au and Ni + Ni reactions were detected 
by the FOPI detector which is a modular detection system with almost Air 
coverage. Of particular importance for the present investigation is the central 
drift chamber CDC which is mounted inside the superconducting magnet with a 
solenoidal field of 0.6 T strength. The CDC has complete cylindrical symmetry, 
it covers the polar angles d from 32° to 150° in the laboratory frame and allows 
the particles masses and charges to be simultaneously determined. For charged 
pions the geometrical boundaries of the CDC guarantees a detection efficiency 
of better than 65%, in particular the midrapidity region \Y^'\ < 0.1 is fully 
covered above a transverse momentum limit p^ > 0.65 (for the definition of the 
variables used here, see below). Notice, however, that to distinguish the tt + from 
protons an additional momentum limit p < 0.65 GeV/c was introduced This 
limit which corresponds to p[ ^ = 3.2 at midrapidity, has no consequences for 
the situation shown in Fig.l, and it is not required for the ir~ . For protons the 
geometrical CDC boundaries 32° < $ < 150° reduce their detection efficiency to 
approximately 45%. 

The situation for particles produced in symmetric reactions at 1.93 AGeV 
bombarding energy is depicted in Fig.l. The fully drawn curves display the lower 
and upper geometrical limits i?i = 32° and $ u = 150° in the laboratory frame, the 
dotted curves are drawn for pions respectively protons with cm kinetic energies 
of 300 MeV. The shaded area displays the momentum distribution of identified 
A(1232) resonances, i.e. the decay proton and pion were measured by the CDC. 
The resonances are Monte Carlo generated using a fixed mass m&. = 1232 MeV/c 2 
and an isotropic Boltzmann distribution with temperature T = 173 MeV, which 
corresponds to a mean kinetic energy of 300 MeV. Because of the decay into 
a correlated (p, 7r) pair the A(1232) resonance can be observed in regions of the 
momentum space which are prohibited for single pions or protons by the $i = 32° 
cut. At midrapidity the A(1232) resonance is detected for transverse momenta 
p{ > 0.5. Notice that the symmetry law 

/(6,$) = /( 7r -e,7r + $) , (1) 

which holds for all symmetric reactions guarantees that the region > can 
be reconstructed by reflection around midrapidity from the region < 0. In 



4 




Figure 1: Left panel: The reconstructed momentum distribution of thermal 
A(1232) resonances with a temperature of 173 MeV. The full curves display 
the geometrical limits of the CDC, the dashed curves display the momenta of 
pions, respectively protons with a kinetic energy of 300 MeV. Right panel: The 
efficiency of reconstructing resonances with masses m-Q and with kinetic energies 
E kin from (p, 7r) pairs detected by the CDC 



this way the properties of around 50% af all A(1232) resonances are known. The 
remaining ones which cannot be reconstructed decay into a (p, tc) pair, where 
either one of the two particles or both are not detected within the geometrical 
boundaries of the CDC. Under the assumptions made, i.e. isotropic resonance 
distribution and isotropic resonance decay, the efficiency of the CDC for detect- 
ing resonances can be calculated as a function of the mass and the kinetic energy 
of the resonance. The efficiencies to detect resonances with masses m-Q and ki- 
netic energies are plotted in the right panel of Fig.l, without employing 
symmetrization. The detection efficiencies decrease rapidly with mass for small 
kinetic energies. Without the proper efficiency corrections, which in this case are 
solely due to the incomplete 4n acceptance of the detector, the contributions from 
small mass values m-Q would be reduced. Notice that from the results published 
in [22] for the Au + Au reaction at 1 AGeV bombarding energy, the expected 
mean kinetic energy of the A(1232) resonance with mass = 1232 MeV/c 2 
is < i?km >= 200 MeV. To obtain the efficiency correction for a given mass the 
data shown in Fig.l have to be weighted by the energy distribution of that mass 
and integrated. 

In addition to the CDC, the information from the forward plastic scintillation 
wall PLA mainly served trigger purposes. The PLA is a high-granularity com- 
ponent of the FOPI detector and covers the polar angles d from 7° to 30° with 
full cylindrical symmetry. It separates particles only according to their charges. 
In the present analysis the particle multiplicities ripLA measured with the PLA 
were required to be larger than 80% of the maximum multiplicity which is defined 
by the lower boundary of the PM5 multiplicity bin, see jl], §]. This multiplic- 
ity criterion selects near-central collisions in the impact parameter range from 
to 0.3 & max , where 6 max is the maximum impact parameter of a given reaction. 
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On the average the fireball then contains around 3/4 of the maximum number 
A = A pvo] + A tSLlg of nucleons. 

In the present study we shall use the identical conventions which were used 
in U], ||. All quantities which refer to the laboratory (lab) or target frames will be 
labeled by small letters. The center of mass(cm) or fireball frames are generally 
characterized by capital letters. Exceptions are those quantities, like the trans- 
verse momentum, which remain unchanged under the transformation of frames. 
For the transverse momentum we sometimes use the quantity p^ where the in- 
dex (°) indicates a normalization by the factor (A ■ P pro j/A pro j) _1 and where A 
is the particle's mass number. Similarly the rapidity Y in the cm frame is nor- 
malized by the factor (Fp ro j) _1 and this quantity is labeled Y^°\ In the present 
experiments the normalization factors for the pion transverse momenta have the 
values 9.47,8.09,7.02 (GeV/cT 1 for 1.06,1.45,1.93 AGeV bombarding energy, 
to obtain the equivalent values for protons these numbers should be multiplied 
with a factor 0.149. For the rapidities the corresponding normalization factors 
are 1.44,1.26,1.12. 

The data from the drift chamber CDC were analyzed by using three different 
tracking algorithms as described in Q. Differences in tracking are mainly ob- 
served in case of the Au + Au reaction, they consist in a slightly different (10%) 
efficiency with which negative pions are identified. Equal tracking efficiencies can 
be obtained by adjusting the limits in which the tracking parameters are allowed 
to vary [|lj. We have verified that our results do not depend on these limits, 
and agree for the different versions of the tracking algorithm. Our final data, 
i.e. transverse momentum and invariant mass spectra, are the average of three 
independent analyses with these different trackers, and this data sample has been 
used before in our study of the pion production in Au + Au U and Ni + Ni [0 
reactions. The accuracy with which the particle momenta can be measured with 
the CDC was determined in |]]|. For the transverse momentum the relative stan- 
dard deviation cr(p t )/p t is 0.04 for p t < 0.5 GeV/c, and then increases to 0.1 for 
transverse momenta p t = 1.5 GeV/c. The accuracy in determining the azimuthal 
angle (p is better than 1°, the accuracy to determine the polar angle $ is better 
for protons (cr($) = 3°) than for pions (cr($) = 5°). These values may be used to 
estimate the accuracy with which the mass m^ of a baryon resonance B can be 
determined. Without going into the details which will be more fully presented in 
the next section, the accuracy estimate is 



where m is the sum of the proton and pion masses and a is the angle between 
p p and p^. Taking representative values, i.e. 



cr(m B ) 




x 
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^2 = 0.06, ^>=0.04, <x(a)tana = 0.09 

Pp Pn 

one obtains for, e.g. the mass of the A(1232) resonance an error of a(m^) = 
20 MeV/c 2 , with decreasing accuracy for larger masses TOb- This estimate is 
certainly sufficient to determine the resonance mass distribution, it is however 
not representative, since a second source of systematic errors comes from the 
reconstruction of the combinatorial background of uncorrelated (p, 7r ± ) pairs, as 
will be discussed in the next section. In the present analysis the number of directly 
measured (p, n^) pairs is larger than 10 6 for a given reaction, to reconstruct the 
background for this reaction we have used approximately 10 times more pairs. 
The statistical errors are therefore not the decisive factor in our analysis. 



3 Methods of analysis 

We have employed two independent procedures to determine the mass distribu- 
tion of baryon resonances excited in near-central collisions of identical nuclei. In 
the first method the pion p t spectrum is defolded to yield the mass distribution, 
in the second method the invariant mass distribution of the correlated (p, 7r ± ) 
pairs is deduced from the data. Both methods rely on the assumption that the 
resonances involved in the production of pions decay via the emission of only 
one pion. This assumption is not necessarily fulfilled as was pointed out in the 
introduction. If the two-pion decay is present it effects the results of the two 
methods in different ways. In the first case the presence of the two-pion channel 
will cause a distortion of the mass spectrum, in the second case it will mainly 
contribute to the background. It is therefore not expected that both methods 
should yield identical results, however these results should be reasonably close 



since the two-pion channel is expected to be weak ||11|| . We therefore also take 
a reasonable agreement as an indication that systematic errors are sufficiently 
small. 

3.1 The analysis of the p t spectra 

The defolding of the experimental p t spectra of charged pions was the earliest 
method applied by Brockmann et al. [T2|, |TJ[ to deduce the mass of the A(1232) 



resonance. Recently this technique was also applied by |23| to defold the p t 
distribution of tt + from Au + Au reactions. It is well known that the p t spectrum 
of 7r + differs from the p t spectrum of ix~ for momenta p t < 0.3 GeV/c, and that 
this difference is caused by the Coulomb interaction between the pions and the 
nuclear matter distribution 0]. In the appendix 6.1 we shall explore this effect 
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Table 1: Temperatures T and flow velocities (3 used in defolding the measured pt 
spectra of tt ± . The r 1 * 1 and values give the measured, respectively calculated 
(eq.(6)) fraction of correlated (p, 7r + ) respectively (p, n~) pairs, is the prob- 
ability that the proton escapes unscattered. The first row is for the Au + Au 
reaction, the next three rows for the Ni + Ni reactions 



energy 


T 


P 




'"e 


k p 


r 


r 

e 


K p 


(AGeV) 


(MeV) 




(%) 


(%) 




(%) 


(%) 




1.06 


81 


0.32 


0.75 ± 0.25 


0.75 


1.0 ± 0.3 


0.6 ± 0.2 


0.22 


2.7 ± 0.9 


1.06 


79 


0.23 


1.0 ± 0.3 


2.17 


0.46 ± 0.15 


0.75 ± 0.25 


0.68 


1.1 ± 0.35 


1.45 


84 


0.29 


0.95 ± 0.3 


2.16 


0.44 ± 0.15 


0.6 ± 0.2 


0.71 


0.85 ± 0.3 


1.93 


92 


0.32 


1.05 ± 0.3 


2.19 


0.48 ± 0.15 


0.6 ± 0.2 


0.71 


0.85 ± 0.3 



more closely and show that the Coulomb effects cancel when the p t spectra of 7r + 
and 7T~ are combined. We follow this prescription in our analysis and add the 
p t spectra of 7r~ and 7r + , measured in a narrow rapidity gap \Y^\ < 0.1. The 
individual p t spectra of it~ and tt + have almost identical cross sections for p t > 0.3 
GeV/c 0, this part of the p t spectra thus supplied the required normalization 
factor. The pion p t spectrum is defolded under the hypothesis that it is entirely 
due to the one-pion decay of excited baryon resonances: 

d<7j " d<7ij(T, (3) 

— = }^ — J - /(mi) ; max j > n , 3 

apt ~i dp t 

where da^/dpt is the midrapidity transverse momentum distribution of pions 
which are emitted in the decay of a baryon resonance with mass m\ > m^. In 
[2~5|1 the method is described to calculate da^/dp t . This distribution depends only 



weakly on the kinetic energy distribution of the resonances, a phenomenon which 



was first noticed by Brockmann et al. [12, 13]. Using the formula of Siemens 



and Rasmussen |24[j the kinetic energy distribution is usually parametrized by 
the nuclear temperature T and the flow velocity (3. For these we have used in 
our analysis the values published in ]22| for the Au + Au, and in [|2£| for the 



Ni + Ni systems. These values of T and (3 are listed in table 1. The defold- 
ing of eq.(3) was performed with the EM algorithm ||26|| . The EM algorithm 



finds after a sufficient number of iterations the Maximum-Likelihood estimate of 
/(mi). Its performance and accuracy, particularly when applied to defolding the 
pion p t spectra, was studied in detail in J27J. It was shown that the result is 
unambiguous and independent of the initial /(mi) distribution, when the func- 
tions cfoy(T, /3)/dp t are linear independent. The mass resolution obtained by the 
defolding was studied as a function of the number of iterations by means of the 
point-spread-functions of m\. Allowing for infinite number of iterations the mass 
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resolution reaches for all m\ its ideal value of 10 MeV/c 2 in accordance with the 
width of the used mass bins. Limiting the number of iterations the mass resolu- 
tion decreases especially for larger masses. The number of iterations is, however, 
bounded by the statistical errors of the measured da-Jdp t spectrum. The it- 
erations are stopped when the calculated spectrum (right hand side of eq.(3)) 
becomes identical with the measured da^/dp t spectrum within these errors. Oth- 
erwise the algorithm would run into describing the statistical fluctuations. The 
mass resolution reached for masses in the important region m ; < 1400 MeV/c 2 
was 10 — 30 MeV/c 2 , for masses m\ < 1700 MeV/c 2 it still remained below 50 
MeV/c 2 . These errors are the statistical uncertainties for a mass bin m u the ac- 
curacy to determine the mean value of the mass distribution /(m;) is much better 
and in fact mainly determined by systematic errors which we have deduced from 
a comparison between our two different techniques to obtain /(mi). 

3.2 The analysis of the (p, 7r ± ) pairs 

The invariant mass of a (p, n ± ) pair can in principle be calculated from the 
measured momenta p p and p^ via 



m B = J (e p + e w ) 2 - (p p + p^) 2 where (4) 



e p = V (Pp) 2 + ( m p) 2 and e vr = V (Pir) 2 + ("O 



\2 



For many of such pairs the resulting mass distribution dn, meas J 'drriB contains 
a large combinatorial background contribution drifo ack J 'dm-Q from uncorrelated 
(p, 7r ± ) pairs, which has to be subtracted from the former in order to obtain the 
required mass distribution dri^/dm-Q of truly correlated pairs: 



= rfw fmeas) _ A _ r± X ^ 

dm-R dm-R ^ ' dm 



_fl_ r ±)_^0. ( 5 ) 



In this expression dn, meas J 'dm^ and dn^^J 'dm-Q are normalized to equal inten- 
sity, and r ± is the ratio of the number of "true" pairs to the number of "random" 
pairs. It is assumed that r ± is constant, i.e. does not depend itself on me, which 
follows from the condition that dri^^J dm-Q is correctly known. The analysis 
consists of determining the size of r^ 1 and the shape of the background contribu- 
tion dnf hack) /dm B . 

The generally accepted way to obtain the background spectrum is based on 
the method of event mixing. This method takes the proton and the pion of a 
pair not from the same event but from different events, using the fact that all 
"true" correlations do not exist in mixed pairs, but that "random" correlations 
will survive the mixing. Besides the true correlations due to the resonance decay 
there are additional correlations due to the reaction dynamics as, for example, 
the particle focussing with respect to the event plane. These latter correlations 
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***** 



Figure 2: The signal for the existence of correlated (p, 7r + ) (left panel), respec- 
tively (p, 7r~) pairs (right panel) in the reaction Au + Au at 1.06 AGeV bom- 
barding energy. The total intensity in the correlation spectra is normalized to 
zero 



are also destroyed by event mixing and therefore are not removed when subtract- 
ing the mixed spectrum from the measured spectrum unless special precautions 
are taken. These require that for the pairs measured in a certain event the 
corresponding mixed pairs come from events with identical proton and pion mul- 
tiplicities and identical orientation of the event plane. The latter requirement 
can always be fulfilled by the proper rotation of both events around the beam 
axis. The condition of equal multiplicities in the measured and mixed spectra 
also ensures that correlations due to pion rescattering in spectator matter [1], 
are properly subtracted. 

The general question of how reliably the uncorrelated background spectrum 
can be reconstructed by means of the event mixing technique was studied in detail 



by L'Hote P9[ . In the majority of cases the difference between the measured and 
mixed spectrum (right hand side of eq.(5)) includes, besides the invariant mass 
distribution from true pairs, a residual contribution from random pairs which 
becomes the smaller, the smaller the differences between the phase space distri- 
butions of protons and baryon resonances. We have confirmed these conclusions 
with our Monte Carlo studies, and have found that in the present case, where 
the decay is dominated by the A(1232) resonance and where the temperatures T 
are moderate, the contribution of random pairs to the right hand side of eq.(5) 
is small. In addition we take the agreement with the results from the defolding 
procedure as independent evidence that the event mixing technique is reliable in 
the present analysis. 

With the method at hand to determine dn^^/drriB it is easy to prove that 
the measured dr& megs -. / dm B spectrum contains a contribution from correlated 
(p, 7r ± ) pairs. For this purpose we assume r ± = and calculate dn^ /dm-Q accord- 
ing to eq.(5). The resulting spectrum has total intensity zero, the correlations 
appear as regions of positive respectively negative intensities for certain values of 
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m B . This is shown in Fig. 2 for the Au + Au reaction at 1.06 AGeV bombarding 
energy. The correlation signals are strikingly strong and they provide the undis- 
putable evidence that part of the (p, 7r ± ) pairs are correlated due to the decay of 
baryon resonances. 

In order to extract the mass distribution of these resonances the ratios r* 1 
have to be experimentally determined. The determination is done in principle, 
i.e. when statistical fluctuations are negligible, by the requirement that the right 
hand side of eq.(5) should be positive. In reality these fluctuations cause part of 
the dn^/drriB spectrum to attain negative values. From Monte Carlo calculations 
with similar counting statistics as in the experiments we have concluded that the 
negative part may vary between 2% and 5% of the total positive intensity. These 
fractions provided the proper criteria for the background subtraction and the 
determination of the r^ 1 values within an uncertainty of 30%. Furthermore, the 
resulting rfn ± /rfm B distribution depends only weakly on the chosen value of 
in the region where dn^/dm-^ has its maximum. The values obtained for r + 
respectively r~ are listed in table 1. 

The relative contribution of correlated pairs r^ 1 can be calculated from the 
measured proton multiplicity n p and the known detection efficiency for protons 

r ± = ^afn% , (6) 

where af is an isospin dependent factor which takes into account that part of the 
detected pions are correlated with undetected neutrons, K p is the probability that 
the proton reaches the detector without another proton - nucleon interaction, and 
£„- describes the fraction of the pions which are correlated, i.e. due to resonance 
decay. The factor af depends on the isospin of the decaying resonance and is 
different for 7r~ and n + . For isospin 1=1/2 one obtains ay 2 = 1 and cxy 2 — 
for 7r~ and 7r + , respectively, for I = 3/2 the corresponding factors £1X6 ^3^2 — 
(N + 2-Z)/(10-N + 2-Z) and a+ /2 = 9 • Z/(10 ■ Z + 2 • N). Notice that I = 1/2 
resonances cannot be observed in the (p, tt + ) channel. The isospin factors for 
1 = 3/2 resonances depend on the N over Z ratio. For the Au + Au reaction 
one obtains a^ 2 = 0.69, a^ 2 = 0.21, for the Ni + Ni reactions these factors are 
a t/2 = 0-74, «3/ 2 = 0.24. Assuming a negligible contribution from thermal pions 

= 1) and that only the decay of / = 3/2 resonances contribute to the yield 
of correlated pairs (af = a0 2 ), eq.(5) provides an upper limit rf of when 
Kp — 1. The values of rf- are quoted in table 1. The direct comparison between 

and rf gives the experimental values of also quoted in table 2. The errors 
of /tp orginate from the uncertainty of (30%). In the next section we shall 
discuss the consequences which follow from these n p values. 



11 



3 

T5 

O 
T3 



o Au 1.06 AGeV 

♦ Ni 1 .06 AGeV 

□ Ni 1.45 AGeV 

A Ni 1 .93 AGeV 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.! 

p, (GeV/c) 



Figure 3: The measured (open symbols) and calculated (full curves) transverse 
momentum spectra in the rapidity range \Y^\ < 0.1 for the reactions Au + Au 
(Au) and Ni + Ni (Ni) at the stated incident energies 

4 Experimental results 

4.1 Results from the p t distributions 

The pion p t spectra obtained by adding the measured 7r + and tt~ spectra in the 
rapidity range —0.1 < < 0.1, are displayed in Fig.3. For the defolding by 
means of the method presented in subsection 3.1 we have used the temperatures 
T and the radial flow velocities (3 determined by [^j and |28] and listed in table 
1. The mass distributions are shown as histograms in Fig.4. The shaded areas 
represent the uncertainties of the defolding technique. These uncertainties are 
firstly due to the possible changes in the values of T and (3, and secondly due 
to the possible contribution of the 2tt decay channel to the low-momentum part 
of Fig.3. Considering the first point first it should be mentioned that a recent 
analysis []30| of the reaction Au + Au at 1.06 AGeV incident energy indicates a 
lower value of the temperature (T « 50 MeV) and a larger value for the collective 
velocity (/3 ~ 0.5) yielding a 10% larger mean kinetic energy of the A(1232) 
resonance than deduced from table 1. In order to estimate the changes which 
may result from different values of the temperature and the flow we have repeated 
the defolding of all pion p t spectra using a 40% smaller value of T and values of 
f3 which correspond to a 10% increase in mean kinetic energy. The uncertainties 
of /(mj) are given by the shaded areas in Fig.4. The uncertainty contributions 
from the decay into two pions were estimated by starting the defolding of the p t 
spectra first at momentum values pt = 0.1 GeV/c and then at p t = 0.15 GeV/c. 
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Figure 4: The invariant mass spectrum of baryon resonances excited by the 
reactions stated in each picture. The shaded areas correspond to the analysis of 
the measured transverse momentum spectra of tt ± , the full points to the analysis 
of the measured (p, tt^) pairs. The arrows point to the maximum of the free 
A(1232) mass distribution 



Also these changes are included in the shaded areas in Fig.4. 

The main features of the mass distributions obtained by the defolding tech- 
nique are their shifts of approximately —60 MeV/c 2 away from the mass of the 
free A(1232) resonance, shown by the arrows in Fig.4, and their weak extensions 
towards larger masses. These tails of the mass distributions are observed for all 
reactions, but it appears to be strongest in case of the Ni + Ni reaction at 1.93 
AGeV bombarding energy. 



4.2 Results from the (p, 7r ± ) pairs 

We first consider the effects which the residual Coulomb interaction between 
charged pions and the nuclear matter distribution might exert onto the recon- 
struction of invariant masses. In Fig. 5 the invariant mass distributions obtained 
from (p, tt + ) respectively (p, ir~) pairs are shown for all reactions. These distribu- 
tions are normalized to equal integrated intensity to allow for an easy comparison 
between the n + and tt~ data. The values of r ± deduced from this analysis and 
the corresponding values of are listed in table 1. In all reactions the values 
of ftp £1X6 larger than those of k+ . This is the behaviour one would expect if the 
resonance mass distribution contains an I = 1/2 component with a mean free 
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Figure 5: The invariant mass spectra of correlated (p, vr ) (open points) and 
(p, 7r + ) (full points) pairs for the reactions stated in each picture 



mass above the A(1232) mass. The large value of n~ in case of the Au + Au 
reaction, c.f. table 1, requires a 25% contribution of / = 1/2 resonances to reduce 
it to its regular size of « 1. In case of the Ni + Ni reactions the k„ 

are of size ~ 0.6 with a 20% contribution of / = 1/2 resonances. These apparent 
differences between the Au + Au and Ni + Ni reactions are not understood yet. 
On the other hand, in all reactions the maxima of the mass distributions shown 
in Fig.5 are shifted towards masses below the A(1232) mass. In order to quan- 
tify these results without taking resort to the Lorentz shape of a free resonance, 
we have calculated the mean mass m-A =< > and the standard deviation 



a a = y < Trig > — < ttib > 2 in the region of the prominent peak with an upper 
limit m-Q < m < . For m < we have chosen a value m< = 1300 MeV/c 2 which is 
motivated by a recent theoretical calculation of the A(1232) mass shape ||3f| and 
which shall be further discussed in the next section, c.f. Fig. 8. The results for 
and <ta are quoted in the first two rows of each reaction entry in table 2. One 
notices that from (p, ir + ) pairs is on the average slightly larger than from 
(p, tt~) pairs, and that the opposite holds for a&. We interpret the first result as 
caused by the differences in the Coulomb interactions between ir + respectively 
7r~ and the baryons, whereas the second result again indicates the existence of 
1 = 1/2 contributions which can only be present in the (p, 7t~) channel. 

The influence of the Coulomb interaction is studied in more detail in the 
appendix 6.1, where it is shown that similar to the prescription employed in case 
of the defolding method, also in case of the (p, ir ) correlations the influence can 
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be minimized by adding the invariant mass spectra of (p, n + ) and (p, ir~) pairs 
to yield the mass spectrum of (p, 7r ± ) pairs. On the latter spectrum we will base 
our further analysis. 

The mass spectra from (p, n ± ) pairs are shown, together with the correspond- 
ing spectra from the defolding method, as points in Fig.4. For each reaction 
these mass spectra obtained by two different techniques are very similar in the 
range around the prominent A(1232) peak. Notice that the efficiency corrections 
represented by Fig.l were not applied to the (p, tt^) data, but other corrections 
discussed in the appendix 6.2 were applied. The reason that the influence of the 
detector geometry onto the invariant mass does not need to be considered is due 
to our finding from Monte Carlo simulations, that the corresponding decrease in 
the acceptance of small masses is counterbalanced by an enhancement caused by 
the finite detector resolution, see the appendix. Both effects are of similar size 
and therefore cancel within the achieved experimental accuracies. 

4.3 The A (1232) mass shift and higher-lying resonances 

For a detailed analysis of the results presented in Fig.4 we have applied the same 
methods used in the analysis of Fig. 5. The mean mass m& and the standard 
deviation a a were deduced in the region m < me < m<, their values are quoted 
in the last two rows of each reaction entry in table 2. The mass shift of Atob = 
—71 ± 6 MeV/c 2 found for the Au + Au reaction appears to be significantly 
larger than the shifts found for the Ni + Ni reactions, which do not seem to 
depend on incident energy and which amount to Atob = —49 ± 4 MeV/c 2 . The 
interpretation of these shifts is postponed to section 5. 

Finally the Fig.4 illustrates that the measured mass distributions extend with 
weak tails to masses above itl-q > m < . Only the (p, tt^) data prove that these tails 
are associated with the excitation of resonances above the A(1232) resonance. 
The defolding technique only proves that this conclusion is not in contradiction 
with the measured pion p t spectra which could also contain thermal pions. The 
proof for the excitation of higher resonances requires that pions and protons 
from the decay of these resonances are correlated. In order to confirm that pions 
with large momenta are indeed emitted by these high resonances and that these 
pions are correlated with protons we have reanalyzed the (p, tt^) data with a 
lower threshold of p t > 0.3 GeV/c for the pion momentum. The invariant mass 
spectra, separated according to the pion charges, are shown in Fig. 6 for the Ni + 
Ni reaction at 1.93 AGeV incident energy. As shaded areas the Fig. 6 also displays 
the invariant mass spectra obtained without the threshold on the pion momenta. 
Two conclusions can be derived from Fig.6: 

• The tail of the mass distribution is strongly associated with high p t pions, 
since the applied p t cut reduces only the yield from the mass region itlb < 
m < . This is, of course, required by the decay kinematics of the resonances, 
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Figure 6: The shaded areas show a logarithmic plot of the invariant mass spectra 
of correlated (p, ir + ) (left panel) and (p, 7r~) pairs (right panel) from the reaction 
Ni + Ni at 1.93 AGeV bombarding energy. The points display the equivalent 
spectra selecting pions with p t > 0.3 GeV/c 

provided such resonances were excited. 

• The shape of the deduced mass distribution is different for (p, 7r + ) and 
(p, 7T _ ) pairs, and the number of correlated pairs has increased more for 
the 7T~ than the tt + . This asymmetry again confirms our earlier conclusions 
about the contribution of I = 1/2 resonances, since these cannot decay into 
(p, 7r + ) pairs. 



5 Summary and discussion 

In the present work we have shown that pions are produced in relativistic heavy 
ion reactions by the decay of excited baryon resonances, with the dominant con- 
tribution from the A(1232) resonance. The heavy ion reactions studied were Au 
+ Au at 1.06 AGeV and Ni + Ni at 1.06, 1.45 and 1.93 AGeV incident ener- 
gies. The mass distributions of the baryon resonances were deduced by means 
of two independent methods: Firstly by defolding the pion p t spectra measured 
at midrapidity by means of eq.(3), and secondly by calculating the invariant 
masses of correlated (p, 7r + ) respectively (p, tt~) pairs by means of eq.(5). In both 
analyses the Coulomb distortions due to the nuclear matter distribution were 
compensated by adding the data obtained for tc + and tt~ . In case of the (p, 7r ± ) 
pairs this compensation was performed with the invariant mass spectra which 
were obtained separately for 7r + and 7r~. The comparison between both suggests 
the additional excitation of I = 1/2 resonances, and it proves the persistence of 
the residual Coulomb distortion in the (p, ir~) respectively (p, n + ) correlations. 

The results from both methods after the removal of the Coulomb distortions 
were shown in Fig. 4 and table 2. The agreement between the mass distributions 
from the two independent methods is the prime indicator that systematic errors 
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Table 2: The first 2 rows of each reaction entry give the values of ira and a a 
deduced from (p, tt + ) respectively (p, n~) pairs, the third row gives the same 
for the average (p, tc^). The last row represents the results of the p t analysis 
including, in addition to m A and <ta, also the fraction tb>a of resonances with 
masses above the A(1232) resonance, and the relative fraction n A of baryons 
excited to the A(1232) resonance 



reaction 


m A 




?B>A 


n A 




(MeV/c 2 ) 


(MeV/c 2 ) 




(%) 


Au + Au 


1160 ± 10 


38 ± 5 






at 


1149 ± 10 


48 ± 5 






1.06 AGeV 


1154 ± 10 


43 ± 5 








1154 ± 5 


48 ± 3 


0.08 ± 0.03 


6.5 ± 2 


Ni + Ni 


1173 ± 10 


41 ± 5 






at 


1171 ± 10 


53 ± 5 






1.06 AGeV 


1173 ± 10 


47 ± 5 








1175 ± 5 


50 ± 3 


0.10 ± 0.03 


11.5 ± 3 


Ni + Ni 


1177 ± 10 


41 ± 5 






at 


1185 ± 10 


56 ± 5 






1.45 AGeV 


1181 ± 10 


48 ± 5 








1176 ± 5 


51 ± 3 


0.15 ± 0.05 


16 ± 4 


Ni + Ni 


1174 ± 10 


40 ± 5 






at 


1166 ± 10 


56 ± 5 






1.93 AGeV 


1171 ± 10 


48 ± 5 








1182 ± 5 


49 ± 3 


0.21 ± 0.10 


22 ± 6 



in our analysis are small. Relative to the mass of the free A(1232) resonance the 
maxima of the deduced mass distributions are shifted towards smaller masses. 
The mass shift is on the average Am A = —60 ± 10 MeV/c 2 , where the lower 
limit corresponds to Au + Au, the upper limit to Ni + Ni near-central reactions. 
The width of the A(1232) remains almost unchanged. The mass shift is more 
than 2 times larger than found in the preliminary analysis of ||20|| . The main 
reason for the deviation is that (p, n ) pairs with small relative angles between 
the particle momenta are not excluded in the present analysis as they were in 
20fl , see also the discussion in the appendix 6.2. Our analysis of the (p, 7r ± ) pairs 



suggests that at least 50% of the detected pions are correlated with protons, i.e. 
they originate from the decay of baryon resonances. For the remaining part, 
the correlations are either destroyed by the residual scattering of the proton by 
nucleons, or those pions are of direct origin. The agreement between the mass 
spectra from the defolding and correlation analyses makes the former possibility 
more likely. In particular the application of p t cuts in the correlation analysis 
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indicates that firstly higher resonances then the A(1232) resonance are involved 
in the pion production, and that secondly at least part of these resonances have 
isospin / = 1/2. 

The mass shift deduced presently and the almost unchanged width of the 
mass distribution are in close agreement with the recent results determined from 
(p, vr + ) pairs in central Ni + Cu collisions at 1.97 AGeV bombarding energy in 



19|| . Nevertheless it should be pointed out that the resonance formula of |3T 



used in the analysis of [19] to describe the A(1232) mass distribution produces a 
tail which overestimates our measured p% spectrum as demonstrated in the upper 
part of Fig. 7. Notice that the inclusion of even higher- lying resonances would 
only enhance the disagreement at large transverse momenta. We thus conclude 
that under the present conditions the mass distribution assigned to the dominant 
A(1232) resonance cannot be adequately described by the parametrization of 



There are two conceivable causes for these modifications of the resonance 
shape: The resonance masses are shifted because of their nuclear environment 
or /and the resonances are in thermal equilibrium with the hadronic matter at a 
low temperature. 

The effect of the hadronic environment on the masses of hadrons was theoreti- 
cally studied in many publications, for a recent review c.f. [02 . The environment 



causes in general a mass shift which can be either positive or negative and depends 
on the hadronic density. The verification of this prediction is at the moment of 
considerable interest, at GSI energies the experiments at present mainly focus 



on the kaon mass j33], [34]]. The situation with respect to baryons appears less 
understood. In a recent study the mass of the A(1232) resonance was calcu- 
lated to be shifted by Am^ m — 10 MeV/c 2 when corrections due to the nucleon 
interaction in the 7rN loop of the A self energy are taken into account. The weight 
function /a(wb) which replaces in this calculation the free mass distribution of 
the A(1232) resonance is used in our further analysis. 

On the other hand a hadronic environment at low temperature T and in 
thermal equilibrium with the baryon resonances will result in an apparent mass 
shift when these resonances decay. Of course, the hypothesis that all hadrons 
are in thermal equilibrium needs verification. But even if equilibrium is not 
established the development of nuclear radial flow will cause the resonances to be 
populated predominantly at the low- mass side in late stages of the reaction ||. 
In the following we will assume for reasons of simplicity that the apparent mass 
shift of the A(1232) resonance can be described by a thermal weight function of 
Boltzmann type with temperature An estimate of Ta, which henceforth shall 
be called temperature although it should be considered primarily as fit parameter, 
is obtained by fitting the mass spectra in Fig.4 by /a(^b) • exp(— (me — wio)/2~a) 
in the mass region m-Q < m < , where /a(^b) was taken from the calculations 
of f35fl. The mass distributions obtained by both methods enter the fits with 



equal weights. The mass range 1150 MeV/c 2 < me < 1250 MeV/c 2 is however 
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Figure 7: Upper part: The p t spectrum of ^ measured(circles) and calcu- 
lated(dotted curves) using the /a(^b) distribution of [19| and the T, (3 values 
published in |[28|| . Lower part: The decomposition of the p t spectrum into the 
A(1232) component (mg < 1.3) and the higher resonances component (mg > 1.3) 



considered with a larger statistical weight, since the very low masses me < 1150 
MeV/c 2 are expected to be stronger affected by the possible systematic errors of 
both analysis methods, and the larger masses me > 1250 MeV/c 2 are stronger 
affected by the contribution from higher resonances. As an example the left panel 
of Fig. 8 displays the fit to the invariant mass distribution in case of the Ni + Ni 
reaction at 1.93 AGeV bombarding energy. The required temperature is in this 
case Ta = 59 ± 8 MeV, for the other reactions the temperatures are listed in 
table 3. The errors represent the uncertainties of the fits deduced by slightly 
shifting the fit region. Notice that the deduced values of Ta depend on the used 
parametrization of the mass distribution. The parametrization of |31| leads to 
values for Ta smaller by approximately 30% than the ones listed in table 3. 

The temperatures Ta obtained in this way are of the order of 50 MeV and 
therefore much smaller than the temperatures deduced for the Au + Au and 

| and listed in table 1. As pointed 
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Ni + Ni reactions by other means 
out earlier in JIJ, these small values of Ta are linked to the enhancement of the 
pion spectrum at very small transverse momenta, which cannot be explained by 
the mass distribution of the free A(1232) resonance, but requires the effective 
mass to be strongly shifted towards the threshold mass by a medium with low 
temperature. 

We interprete the calculated mass distribution of Fig. 8 as the sole contribution 
from the A(1232) resonance, which for m-Q = m < = 1300 MeV/c 2 practically has 
reached zero intensity. The remaining yield not accounted for is, because of 
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Table 3: Temperatures Ta and radial flow velocities (3& deduced from a fit to 
the mass distributions in the range m-Q < 1300 MeV/c 2 , and a fit to the kinetic 
energy distributions, c.f. Fig. 8. The first row is for the Au + Au reaction, the 
next three rows for the Ni + Ni reactions 



energy (AGeV) 


T A (MeV) 


Pa 


1.06 


41 ± 5 


0.46 ± 0.03 


1.06 


53 ± 8 


0.39 ± 0.05 


1.45 


60 ± 8 


0.47 ± 0.04 


1.93 


59 ± 8 


0.53 ± 0.04 



its correlation with protons, at least partly caused by higher lying resonances. 
As pointed out before this interpretation is also suggested by the fact that the 
(p,7r~) and (p, 7r + ) invariant mass spectra are different for large pion momenta 
and that both pairs have different values, see table 1. The upper estimate 
of the ratio between the contribution from high-lying resonances and the one 
from the A(1232) resonance is quoted in table 2 as tb>a- These values were 
extracted from the results of the defolding technique, since the analysis of the 
correlation results, also shown in Fig.4, requires information about the admixtures 
of / = 1/2 and 1 = 3/2 states in the tails of the invariant mass distributions 
which does not exist with the required accuracy. We have used the latter results, 
however, to obtain the quoted errors of tb>a- The contribution of the higher 
lying resonances to the p t spectrum, obtained in this way, is shown in the lower 
part of Fig. 7. It demonstrates the difference to the alternative where the high- 
momentum tail is interpretated as mainly due to thermal pions: The distribution 
from thermal pions would continue to rise towards small p t values and would give 
a thermal contribution which is larger than the values of tb>a quoted in table 2, 



in agreement with [35 



With these informations and the known [l], |2| dependence of the pion mul- 
tiplicity on the number of participants A part , one may calculate the fraction ha 
of participants which are excited to A(1232) resonances under the present ex- 
perimental conditions. The values of ba are listed in table 2, they increase with 
bombarding energy primarily because the number of pions increases. Similar val- 
ues for the fraction of A(1232) resonances populated in the Ni + Ni reactions 
were extracted in ||28j| , where the high-momentum tail was interpreted as due to 
thermal pions emitted from an equilibrated radially expanding source with the 
temperatures T and the average expansion velocities f3 listed in table 1. 

According to the observed mass shifts in the current analysis the freeze out 
conditions for pions from late A(1232) decays would be characterized by a much 
lower value of Ta = 40 — 60 MeV. This value when interpreted as temperature is 
of similar size as the average freeze out temperatures of protons in the complete 
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Figure 8: Left panel: The measured (symbols) invariant mass spectrum of the 
Ni + Ni reaction at 1.93 AGeV bombarding energy, and the fit (curve) using 



the Boltzmann corrected /a(^b) distribution of | 3q| . Right panel: The kinetic 
energy spectrum (symbols) of the A (1232) resonances from the same reaction. 
The curve is a fit employing the Siemens- Rasmussen formula |24| . The values for 
the fitted temperature (left panel) and radial velocity (right panel) are listed in 
table 3 



phase space, extracted in recent analyses of the same present data samples |30], ^6 
which also point to a larger collective energy and thus smaller temperatures than 
estimated from the kinetic energy spectra of pions and light charged particles 
alone. It is also in line with the analysis of the Au + Au reaction at lower 
bombarding energies 



On the other hand, the assumption of thermal equilibrium and instant freeze 
out at Ta fails to reproduce the A(1232) yield. As was shown in P3] , where the 
attainment of equilibrium was postulated, the considerably higher temperatures 
listed in table 1 are necessary to describe the observed production probabilities. 
These findings are in agreement with [331 and are confirmed by an improved 



calculation where the widths of all higher-lying resonances were included via a 
Breit-Wigner parametrisation with parameters listed in the particle data booklet 
[39Q . For a consistent treatment of the yields and the energy distributions of 
baryons and pions we drop the surface corrections and use the phase space density 
of infinite nuclear matter |H| . At a temperature of T = 92 MeV this model then 



predicts a value of approximately tb>a = 0.35 for the ratio between higher-lying 
and the A(1232) resonance, in variance with the measured value given in table 2. 
From the thermal model one may also deduce the contribution of thermal pions 
which increases with decreasing freeze out density. Estimating the upper limit 
of this contribution from the value of table 1 and the published pion yield 
[[[]. the freeze out density p should be close to normal nuclear matter density p , 
whereas it is predicted that the majority of pions, i.e. those with small momenta, 
freeze out at p « p Q /3 (40], |Ij. 



These inconsistencies in the interpretation of the different experimental ob- 
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Figure 9: The left panel displays the calculated p t distributions of it , n~ , and 
7T + (curves). The circles present the sum of the 7r~ and 7r + distributions. The 
centre panel displays the invariant mass distributions obtained by defolding the 
p t spectra of the left panel. The right panel illustrates the reduction of low- 
momentum 7T° caused by their absorption in the source 



servations probably point to the fact that the simple picture of an instantaneous 
freeze out in thermal and chemical equilibrium is not justified. The number of 
baryon resonances and pions is determined very early in the reaction when the 



energy to excite these resonances is still sufficiently large ||10|| . The dominant 
part of the A(1232) resonances is measured at late stages of the reaction, when 
this energy has considerably decreased because of earlier it emissions, and be- 
cause of the development of radial flow. In case one is allowed to interpret Ta 
as that part of the energy corresponding to unordered motion at freeze out one 
may deduce from resonances in the mass range me < m < also the amount of en- 
ergy converted into flow at the late stage of the reaction. We have reconstructed 
for those resonances their kinetic energy spectrum which is depicted in the right 
panel of Fig.8 for the Ni + Ni reaction at 1.93 AGeV bombarding energy. The 
reconstruction procedure is equivalent to the one used to identify the resonances, 
i.e. the E — tub spectrum of mixed pairs is subtracted from the one of measured 
pairs with a ratio r which was determined earlier and is listed in table 1. The 
kinetic energy spectra are corrected for the effects of the CDC geometry onto 
E and m&. The resulting E^ [n distribution of the shifted A(1232) resonance ap- 
pears to have a maximum close to E kin = 0.2 GeV. The shape of this distribution 
suggests the existence of a strong collective flow. The fit by means of the Siemens 
- Rasmussen formula [23 with fixed mass rriA = 1180 MeV/c 2 and temperature 
Ta = 59 MeV and shown as full curve in the right panel of Fig.8, yields a flow 
velocity /3a = 0.53 ± 0.04. The values of /3a deduced similarly from the other 
kinetic energy spectra of low-lying resonances are listed in table 3. If interpreted 
in this way it shows that at late freeze out times of the pions the temperatures 
do not change significantly in the Ni + Ni systems, but the radial flow velocities 
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increase with energy. The Au + Au system, compared with the Ni + Ni system, 
is characterized by a lower temperature and larger radial flow velocity. 

Because of the very low nuclear temperatures deduced from our analysis one 
cannot neglect the possibility that the shifts of the resonance mass might be 
induced by still other effects of the nuclear environment than those already con- 
sidered in P5 |. In |32| the alternatives are discussed in some detail. In case 
the mass shift scales with the nuclear density and higher lying resonances are 
excited more abundantly in the earlier times of large densities the observed mass 
spectrum of the baryon resonances would be rather complex and difficult to de- 
compose in contributions from specific resonances. From Fig. 6 it appears that 
whatever the true mechanism is, it produces a resonance mass distribution with 
near-exponential decline towards larger mass values. It was indeed this expo- 



nential decline which was used very early by Chapline et al. [f|2| to study the 
properties of dense hadronic matter. 



6 Appendix 

6.1 Effects due to the Coulomb interaction 

The Coulomb residual interaction between the nuclear matter distribution and 
charged mesons modifies the primordial phase space distributions of the latter 
after they were emitted by resonance decays. The size of this modification is 
particularly strong for pions because of their small mass, and it is readily observed 
in experiments: The momentum distribution of positive and negative pions differ 
for small momenta, the difference is still observed in systems with total charge 



Z < 60 (pf. In a recent theoretical paper based on BUU calculations |40|| , it was 
proposed that the n^/n + difference provides a sensitive method to investigate 
the nuclear matter distribution at SIS energies. 

In our Monte Carlo studies of how the Coulomb effects may change the de- 
duced mass distributions of baryon resonances we have employed a slightly sim- 
pler model which is based on a spherical and radially expanding baryon distribu- 
tion in thermal equilibrium. The A (1232) resonances with the mass distribution 
/a(^b) are present in the baryon distribution with fraction %, they can decay 
via 7r + , 7r~, or tt° emission only at the surface of this distribution. After the decay 
the emitted particles are propagated in the Coulomb field of the baryon distri- 
bution. The positions and the momenta are calculated in relativistic kinematics, 
the time steps are sufficiently small to handle the complete momentum range of 
pions. The final momentum is reached, when the relative change of the momen- 
tum becomes smaller than 10~ 6 . In case the total energy approaches the pion rest 
mass, this is only possible for 7r~ because of their negative Coulomb potential, 
the particle is considered to be in a bound state and cannot be detected. Further- 
more we use an absorption probability of 0.2 for those pions, which revert into 
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Figure 10: The influences of finite detector acceptance and resolution on the 
resonance mass (left panel) and kinetic energy (right panel) distributions deduced 
from correlated (p, tt°) pairs 



the baryon distribution. Notice that absorption can occur for all pion charges, 
the probability depends on the relation between the momentum of the decaying 
resonance and the momentum of the pion emitted by the decay. To be close to 
the situation of the Au + Au reaction we consider a source of normal nuclear 
matter density with 310 baryons and a total charge given by the N/Z ratio of Au 
[fi"Tll . The kinetic energy of the baryons is assumed to be distributed according to 
the formula of Siemens and Rasmussen |M] with a temperature T = 81 MeV and 
an expansion velocity (3 = 0.32. In addition to the geometrical and the detection 
limits of the CDC we also include the finite momentum resolutions of the CDC, 
characterized by the parameters cr(p t )/p t = 0.05, cr(y?) = 1°, and ct(t9) = 5° 0. 

This simple model is sufficient to reproduce the gross features of the mea- 
sured |l|, |2| transverse momentum spectra of 7r ± , i.e. their difference at momenta 
Pt < 0.3 GeV/c, and their identical decline for higher transverse momenta, see 
Fig. 9 left panel. It is also evident that the sum of the tc + and tc~ momentum 
distributions is a good approximation for the primordial momentum distribution 
obtained from tt°. It is therefore not surprising that the primordial A(1232) 
mass distribution /a(^b) is not retrieved from the distorted pion momenta. The 
defolding of the pion p t spectra yields from the n~ respectively ?t + a /a(^b) 
distribution which is shifted to smaller/larger masses, see Fig. 9 centre panel, 
whereas the sum of both p t spectra allows to obtain with good approximation 
the undistorted /a(^b) distribution as expected since the Coulomb effects cancel 
when the p t spectra of tc~ and tc + are added. 

It might be interesting to illustrate how the pion absorption in nuclear matter 
changes the pion p t distributions. In Fig. 9 right panel the p t spectra of 7r° are 
compared with and without absorption. It is obvious that absorption mainly 
suppresses pions with small transverse momenta. On the other hand this part 
of the momentum distribution becomes enhanced again when the decay channel 
into 2/T opens. Therefore the net effect is presumably small, and it is beyond the 
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Figure 11: Comparison between the resonance mass distributions from (p7r°), 
(j97r~), and (pn + ) pairs. The (p7r ± ) distribution corresponds to the normalized 
sum of the (pn~) and (pn + ) distributions 



scope of the present study to disentangle these effects experimentally. 

The situation becomes more complicated when looking at the correlations 
between pions and protons. In Fig.l we have shown how the geometric boundaries 
of the CDC reduce the detection efficiency for small resonance masses. In the 
left panel of Fig. 10 the uncorrected invariant mass spectrum from (p, 7r°) pairs is 
plotted which shows the size of this effect. On the other hand, when the finite 
detector resolutions are considered one finds an enhancement of the invariant 
mass distribution at small masses. Including both effects one obtains in our 
Monte Carlo simulations an invariant mass distribution very close to the input. 
The cancellation of different effects is not always observed. In case of the energy 
distribution of resonances we find in the right panel of Fig. 10 that these are mainly 
effected by the geometric boundaries, and that corrections have to be applied 
to restore the original distribution from the measurement. The corrections are 
considered reliable only for kinetic energies E^ in > 100 MeV. 

The effects of the residual Coulomb interaction on the invariant masses from 
(j>, 7r ± ) pairs are shown in Fig. 11. Qualitatively the situation is similar to the 
analysis of the pt spectra shown in Fig. 9. The /a(^b) distribution deduced from 
the (p, 7r~) respectively (p, 7r + ) pairs are also shifted to smaller/larger masses 
relative to the mass distribution from correlated (p, tt°) pairs. On the other hand 
the average mass distribution from both (p, tt~) and (p, n + ) pairs comes very 
close to the undistorted distribution. 
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Table 4: The values of the mass m A and the width <7a obtained from tt°, tx , 7r + , 
and the average 7r ± by the analysis of (p, n) pairs or by the defolding technique 





pair analysis 


defolding 




m A 
(MeV/c 2 ) 


(MeV/c 2 ) 


m A 
(MeV/c 2 ) 


(MeV/c 2 ) 




1201 ± 3 


51 ± 2 


1198 ± 2 


48 ± 1 


7T~ 


1197 ± 3 


50 ± 2 


1180 ± 3 


56 ± 2 


7T + 


1207 ± 3 


49 ± 2 


1228 ± 2 


34 ± 2 




1202 ± 2 


50 ± 2 


1202 ± 2 


51 ± 2 



The sizes of the mass shifts m A and the widths <r A deduced from our Monte 
Carlo studies are listed in table 4. It is evident that the results from the corre- 
lation analyses are less affected by the Coulomb distortion than those from the 
defolding analysis. Nevertheless adding the n~ and tt + data one comes closest 
to the primordial conditions, which is the procedure employed in the present 
analysis. 



6.2 Effects due to the tracking efficiency 

The discussion in section 2 is primarily concerned with the efficiency of recon- 
structing the baryon resonances within the geometrical boundaries of the CDC. 
Another efficiency reduction, which requires a careful investigation, is imposed by 
the limited capabilities of the available tracking algorithms to separate and iden- 
tify individual particle tracks in the CDC. These inefficiencies become largest 
when the mass of the decaying resonance is close to the threshold mass mo, i.e. 
when the relative momentum between the proton and the pion is small and the 
two particle tracks are close to each other. For the tracking algorithm, which 
analyzes particle tracks in the xy plane perpendicular to the z direction of the 
magnetic field, the closeness of tracks is given by the relative angle a t = ip^ — ip^ 
and by the magnitude of p[ n \ The smaller a t and the larger p{ the more difficult 
it becomes to separate pion tracks from proton tracks. 

The size of the efficiency reduction as function of a t and p[^ can be obtained 
from a comparison of the invariant mass distributions dn (j neas '> an d dn ^ lxed '> as 

r elms diriB 

functions of a t and p^ since in the latter all pairs are identified with identical 
efficiencies. We show in Fig. 12 the ratio to = ^(mixel) ^ or iP-> 7l ~)-> (p^^)^ an d 
Monte Carlo generated (p, vr) tracks as functions of a t , where the Monte Carlo 
simulation illustrates the ideal case of complete tracking efficiency. In the simula- 
tion the conditions are similar to the experimental ones, in particular the fraction 
of (p, ti ± ) pairs due to resonance decay is in the order of 1%. The comparison 
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Figure 12: Efficiency to identify a proton and a pion track as function of the 
relative angle a t between these two tracks in the x, y plane perpendicular to the 
magnetic field 



with the Monte Carlo simulation proves, that the probability to miss correlated 
(p, vr ± ) pairs never exceeds 20% for a given value of a t , and that it is certainly 
more appropriate to correct for this inefficiency than to remove all data with 
| at I < 50°, as was done in |20| . 

That the efficiency reduction is different for (p,ir~) and (p, tt + ) pairs was to 
be expected since the p and 7r + have the same curvature signs, whereas these are 
different in case of p and 7r~. Furthermore the fact that in case of (p, 7r~) pairs 
the efficiency is reduced only for negative a t angles offers the possibility to check 
the consistency of the efficiency correction. The invariant mass spectra -p— have 
to be the same for a t < and a t > 0. This check was done in our analyses and 
was shown to yield identical distributions within the given uncertainties. 
All data shown in this paper were therefore efficiency corrected according to the 
results of Fig. 12 and the similar dependences on p\ . The sizes of the efficiency 
corrections were directly obtained from the data, their dependence on a t and p\ 
was, however, smoothed by the spline algorithm. 
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